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ABSTRACT 



Dark-matter halos show a universal density profile with a scaling such that less 
massive systems are typically denser. This mass-density relation is well described by a 
proportionality between the characteristic density of halos and the mean cosmic density 
at halo formation time. It has recently been shown that this proportionality could be 
the result of the following simple evolutionary picture. Halos form in major mergers 
with essentially the same, cosmogony-dependent, dimensionless profile, and then grow 
inside-outside, as a consequence of accretion. Here we verify the consistency of this 
picture and show that it predicts the correct zero point of the mass-density relation. 
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1. INTRODUCTION 
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>< 

High resolution cosmological N-body simulations show that the spherically averaged density 
profile of dark-matter halos has the following universal form (Navarro, Frenk, & White 1997, 
hereafter NFW; but Moore et al. 1997) 

p(0 __ Sc m 

with £ the radial distance to the halo center in units of the scale radius r s , and 5 C the characteristic 
halo density p c in units of the critical density for closure p C rit- The parameters r s and p c are linked 
by the condition that the mean density within the virial radius R of a halo of a given mass is a 
constant factor a times the cosmic critical density; here we take a = 200. Therefore, the density 
profiles of halos at a given epoch depend on their mass M through one unique parameter, 5 C or 
x s = r s /R. 

The universality of this profile is not surprising as far as the gravitational clustering of dark- 
matter is self-similar in scale- free cosmogonies and approximately so, over a certain range of masses 
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and times, in other cosmogonies. However, the specific form of this profile and the scaling functions 
5 C (M) or x s (M) at a given time t is hard to predict. Much of the research on this issue has focused 
on the effects of spherical infall (Gunn & Gott 1972; Fillmore &; Goldreich 1984; Bertschinger 1985; 
White <fe Zaritsky 1992; Avila-Reese, Firmani, & Hernandez 1997; Henriksen & Widrow 1998) 
and/or secular evolution owing to dynamical friction and tidal disruption on small mass clumps 
(Syer & White 1997; Nusser & Sheth 1998). Both these mechanisms require a certain stability of 
the system and, therefore, presume the growth of halos governed by minor mergers or accretion. 
The opposite viewpoint that the typical halo density profile is the result of repetitive, asymmetrical, 
major mergers producing strong violent relaxation has also been investigated (e.g., Duncan, Farouki 
& Shapiro 1983). In fact, accretion and major mergers alternate in hierarchical clustering (see, e.g., 
Tormen, Bouchet & White 1997 for the CDM cosmogony). Hence, one should in principle take into 
account the combined action of these two extreme processes. 

NFW have found that the smaller the mass of halos, the denser they are. This mass-density 
correlation was interpreted as reflecting the fact that, in hierarchical clustering, less massive halos 
form earlier when the mean density of the universe is higher. NFW showed, indeed, that the 
characteristic density of halos with a given mass at a given epoch is proportional to the mean cosmic 
density at their time of formation. There is, however, a small caveat in this reasoning. According 
to it the structural properties of halos would be fixed on formation. As major mergers yield a 
substantial rearrangement of the system, while accretion only causes a smooth secular evolution, 
it seems consistent to define the formation of halos as the last major merger they experience, their 
identity being kept during accretion. Yet, the halo formation time used by NFW, drawn from the 
extended Press-Schechter (1974, PS) formalism, does not match this definition. 

Recently, Salvador-Sole, Solanes, & Manrique (1998, SSM) have repeated the same analysis as 
NFW but using a better suited formation time estimate drawn from a modified version of the PS 
theory. This incorporates a phenomenological distinction between accretion and major mergers. 
For a fractional mass capture threshold for major mergers, A m , equal to 0.6, SSM confirm that, 
in any hierarchical cosmogony, there is a clear proportionality between the characteristic density 
of halos and the mean density (or still the critical density) of the universe when they formed. 
Furthermore, SSM have shown that this proportionality is consistent with the following simple 
evolutionary picture for the structure of dark-matter halos. Halos form, as the result of major 
mergers, with essentially universal (i.e., independent of mass and time), cosmogony-dependent, 
scaling parameters S c and x s . Then, until the next major merger, the corresponding dimensional 
parameters p c and r s remain fixed, there only being a shift in S c and x s as the critical density 
decreases and the virial radius expands accordingly. 

The requirement of steady profiles during the accretion phase which follows from the empirical 
scaling and the PS clustering model agrees with accurate dynamical studies of self-similar infall 
(Fillmore Sz Goldreich 1984; Bertschinger 1985; Henriksen, & Widrow 1998). This gives strong 
support to the validity of that evolutionary picture. However, its consistency is still to be checked. 
On the other hand, although this picture might explain the observed proportionality between the 
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mass and characteristic density of halos, the zero point of the mass-density relation in any given 
cosmogony remains to be explained. In the present Letter we address these two issues. 



2. THE PREDICTED STRUCTURE OF HALOS FORMED IN MAJOR 

MERGERS 

In this section we derive, in an arbitrary hierarchical cosmogony, the values of the scaling 
parameters 5 C or x s of halos that result, according to the proposed evolutionary picture, from the 
major merger of their typical progenitors with appropriate structural properties. All elements we 
need for this calculation are given in SSM except for the typical mass (and number) of progenitors 
of a halo with mass Mq at to. This is next derived in the same frame of the modified PS clustering 
model used in SSM. 

The specific capture rate, r c {M' <— M, t), is the rate at which halos of final mass M' capture, 
at t, other halos of mass M < M' . Following Manrique & Salvador-Sole (1996) this is related to 
the specific major merger rate, r m (M —* M',t), given in SSM through 

r c (M' *— M,t) = r m (M -> M', t)^£| , (2) 

where N(M,t) is the mass function. For halos with mass M larger than M'/(A m + 1) the specific 
capture rate is null. This reflects the fact that such halos evolve by accretion and are, therefore, 
identified with the final ones with mass M' which cannot capture themselves. On the other hand, 
the capture of halos with masses M smaller than M'A m /(A m + 1) do not cause, in general, the 
destruction of the capturing halos and, hence, do not imply the formation of new ones. Only those 
captures of halos with masses M ranging between M'A m /(A m + 1) and M'/(A m + 1) necessarily 
correspond to major mergers implying the formation of new halos. Hence, were all major mergers 
binary, the specific capture rate given by equation (pi) would be, within that range, a symmetrical 
function of M around M'/2, which would lead to the equality between the two integrals 

m' m' 
C(M',t)= rf m+1) r c (M' ^M,t)dM= P r c (M' <— M, t) dM. (3) 

T (A m + 1) 

Moreover, the probability that one of the two progenitors of a halo of mass M' has mass M within 
the limits of any of the two preceding integrals would then be given by 

_ r^[M'^MMM',t)] 
9p[M) C(M',t) ' lj 

and the typical masses Mi and M2 of the progenitors of a halo of mass M' would coincide with 
the median value of this distribution function and its complementary value, respectively. N-body 
simulations show that major mergers are binary in agreement with the fact that, in all cosmogonies 
analyzed, the capture rate (Q) is found to be very closely symmetrical around M'/2 in the appro- 
priate range of masses, the fractional difference between the two integrals in equation (||) being 
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less than 3% at any redshift. Thus, the use of the probability (||) to derive the typical mass of 
progenitors is fully justified. 

Let us consider now two gravitationally bound halo progenitors with masses M\ and M 2 at 
turnaround. Each mass Mj includes not only the mass of the main relaxed body, but also that of 
the surrounding matter that will have accreted onto it by the time the merger occurs. Assuming 
radial orbits, the total energy of the system at turnaround is 

£ta = Ui + U 2 + E 12 , (5) 

with Ui the internal energy of each component, and £12 the interaction energy. If we apply energy 
conservation to each progenitor considered in isolation, the internal energy Ui is approximately 
equal to the internal energy of the corresponding fully relaxed halo prior to merger. Assuming 
spherical symmetry and virial equilibrium, this is given by 

1 GMf . . 

with Ri and x s i the virial and dimensionless scale radii, respectively, of progenitor i at merger 
time and F(x s ) a dimensionless function of order unity dependent on the density profile (|l]). The 
interaction energy is also approximately given by the expression 

E* = -™^, (7, 
with D m the turnaround separation between the two centers of mass, equal to 

"t f (M ,to) nl> 



D't = 2GM 



in 



7T 



(8) 



for two point masses falling radially towards each other, where we have taken into account that 
the cosmic time corresponding to zero separation marking the merger of the two progenitors is the 
formation time, if(Mo,to)> of the final halo with mass Mq at to- Finally, taking into account that 
the internal energy U of the final halo at formation, of the form given in equation (|6|), coincides 
with the total energy of the system at turnaround (eq. (5|), we arrive to the approximate relation 

Ffe)+ b£) -j^ 2 - < 9 > 

with M = Mi + M 2 , t equal to tf(Mo,to) i n units of the Hubble time at that epoch, and x s f the 
value of x s at halo formation. 

In deriving the previous relation we have neglected the external gravitational torques acting 
on the system as well as the mutual tidal interaction between the progenitors. External torques 
make the orbits be non-radial which, for a fixed turnaround separation, increases the infall time. 
To follow non-radial orbits it is convenient to use the dimensionless variable S = J 2 /J^(E), where 
J is the specific angular momentum of M 2 relative to Mi, E the corresponding orbital energy, 
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and J£(E) the specific angular momentum of a circular orbit of energy E relative to M\. For an 
unperturbed elliptical orbit of eccentricity e one has 5 = 1 — e. So the values of 5 range from 
(radial orbit) to 1 (circular orbit). Taking the merging time equal to one orbital period leads to 
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*f(M ,t ) 
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(10) 
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instead of equations (0) and (|8|), respectively. This introduces an extra factor in the last term on 
the right-hand of equation (|9|) which depends on the non-null value of 5 and on M 2 /M\. Note 
that, to be more accurate, we should adopt a somewhat increased effective value of 5 to account 
for the fact that the merging time for low eccentricity orbits is typically larger than one single 
orbital period. Concerning the tidal interaction between the two progenitors, this causes the so- 
called tidal braking (part of the interaction force is spent in deforming the systems which slows 
down the orbital motion and increases the infall time) which can also be accounted for by further 
increasing the value of 5. Note that the importance of external torques is expected to increase with 
decreasing mass while, the more massive the halos, the more marked their mutual tides. Therefore, 
both effects can be taken into account by adopting, in a first approximation, a constant positive 
value of 5 in the whole range of masses. 

From the clustering model given in SSM we can compute the values of r, Mi, and M 2 , as well 
as the typical formation time of each progenitor. Then, the evolutionary picture for the structure 
of halos we wish to check allows us to calculate the values of x s i of the two progenitors prior to 
merger from the assumed universal value of this parameter on formation. (This can be done by 
using any of the two independent fitting formulae given in SSM, the two results typically differing 
by less than 10%.) Therefore, equation (||) can be solved for x s f. 



3. COMPARISON WITH THE EMPIRICAL STRUCTURE OF HALOS 

In Figure |l] we plot, for the NFW density law, the fractional difference between the solution 
of equation (|9|), x| f , and the empirical value of x s f obtained by SSM used as the formation value 
for the progenitors, . clS El function of Mo at the present time. Mq is scaled to M* defining a 
unity rms density fluctuation so that similar graphs hold at other epochs. The curves plotted here 
correspond to the standard CDM cosmogony with b = 1.59 although any other cosmogony analyzed 
(in particular those studied in SSM) lead to very similar results. For 5 = 0, there is very good 
agreement between and at the large mass end. The deviation increases towards small Mq 
but, within the four decades mass range shown in this Figure, it stays below 25%. Increasing the 
value of 5 tends to further improve the agreement at small masses without spoiling too much the 
good predictions at the large mass end. 

Moreover, the evolutionary picture analyzed implies one single, cosmogony-dependent, value of 
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x s f fixing the zero point of the mass-density relation. In Table 1 we quote, for various cosmogonies 
(with parameters listed in the first four columns), the value of minimizing, for S = 0.5, the 
maximal fractional difference between x| f and in the four decades mass range analyzed (column 
5), and the empirical value of x^ f obtained in SSM from the x _sc l uare fit of the mass-density 
relation (column 6) with the corresponding 90% confidence interval (column 7). As can be seen, 
the theoretical values favored by the evolutionary picture proposed fall inside the 90% confidence 
interval associated with the empirical values in 6 of 8 cases. Only in the case of the flat, n = —1.5, 
and the open, n = — 1, cosmogonies the theoretical prediction in column (5) is a little too small. In 
fact, the trend for the theoretical value to be slightly smaller than the empirical one seems quite 
general suggesting the existence of some slight bias of the kind mentioned in SSM. In any event, the 
general agreement between the theoretical and empirical values of x s f in all cosmogonies considered 
is hardly a coincidence. This strongly supports the validity of the schematic evolutionary picture 
for dark- matter halos proposed in SSM. 



4. DISCUSSION 

In § 2 we have considered only elliptical orbits with different possible values of S to account 
for their unknown typical eccentricity. These correspond, indeed, to the most reasonable initial 
condition for mergers regardless the cosmogony. Mergers from parabolic or hyperbolic orbits are 
certainly much more improbable because of their non-negative orbital energy. This seems to con- 
tradict the idea that in scale free cosmogonies (i.e., an Einstein-de Sitter universe with power law 
spectrum: P(k) oc k n ) with n = 1, mergers occur typically between halos describing parabolic or- 
bits, while for n > 1 (n < 1) they would follow hyperbolic (elliptical) orbits. The reason is that the 
mean cosmic density and, hence, the mean density of halos with characteristic mass M*, scales as 
M* ( 3+n )/ 2 _ Therefore, for n = 1 one has M* oc R implying that the specific binding energy of such 
halos is constant in time. One must be very cautious with this kind of qualitative arguments. In 
the previous reasoning, the binding energy refers only to a very specific class of halos, those which 
have the characteristic mass at the given epoch. Furthermore, the typical progenitors of such halos 
will hardly coincide with the class of objects which have characteristic mass at the adequate epoch. 

In this Letter we have checked the consistency, relative to energy conservation, of the scaling 
evolution inferred in SSM from the empirical mass-density relation. We must stress however that 
such a relation is not the simple consequence of energy conservation in mergers. Our results show 
that the mass-density relation depends, in addition, on the assumed evolution of halo structure. 
Any departure from the concrete evolutionary picture proposed (e.g., other values of parameters 
5 C and x s at halo formation or a non-steady behavior of halos during accretion) would have led, in 
general, to a mass-density relation in disagreement with iV-body simulations. 

This work has been supported by the Direction General de Investigation Cienti'fica y Tecnica 
under contract PB96-0173. 
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Fig. 1. — Fractional difference between the value on formation of the scaling parameter x s ob- 
tained for halos of mass Mq and the one assumed for their corresponding progenitors (x| f and x| f , 
respectively), for the case of the NFW universal density law and the standard CDM cosmogony. 
The values of the orbital parameter S used range from zero (upper curve) to one (lower curve) in 
steps of 0.25. Solid lines correspond to equal to the empirical value of x s found in SSM. In 
dashed line, the best theoretical prediction obtained (for any fixed value of S) by taking 
free parameter. 
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Table 1. Scaling of Halo Density Profiles 
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